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Constructing  a  Unitary  Hessenberg  Matrix  from  Spectral  Data 

Gregory  Ammar1,  William  Gragg2,  Lothar  Reichel3 

In  memory  of  Peter  Henrici 


Abstract 

We  consider  the  numerical  construction  of  a  unitary  Hessenberg  matrix  from  spectral  data  using 
an  inverse  QR  algorithm.  Any  unitary  upper  Hessenberg  matrix  //  with  nonnegative  subdiagonal 
elements  can  be  represented  by  2n  —  1  real  parameters.  This  representation,  which  we  refer  to  as 
the  Schur  parameterization  of  //,  facilitates  the  development  of  efficient  algorithms  for  this  class  of 
matrices.  We  show  that  a.  unitary  upper  Hessenberg  matrix  //  with  positive  subdiagonal  elements 
is  determined  by  its  eigenvalues  and  the  eigenvalues  of  a  rank-one  unitary  perturbation  of  II.  The 
eigenvalues  of  the  perturbation  strictly  interlace  the  eigenvalues  of  II  on  the  unit  circle. 

AMS(MOS)  Subject  Classification:  15A18,  65F15. 

Keywords:  inverse  eigenvalue  problem,  unitary  matrix,  orthogonal  polynomial. 


Presented  at  the  NATO  Advanced  Study  Institute  on  Numerical  Linear  Algebra,  Digital  Signal 
Processing,  and  Parallel  Algorithms,  Leuven,  Belgium,  August,  1988.  This  research  was  supported 
in  part  by  the  National  Science  Foundation  under  grant  DMS-8704196,  and  by  the   Naval 
Postgraduate    School   Research   Council 


Northern  Illinois  University,  Department  of  Mathematical  Sciences,  DeKalb,  If  60115,  USA 
Naval  Postgraduate  School,  Department,  of  Mathematics,  Monterey,  CA  930  1.'!.  USA,  on  leave  from 
University  of  Kentucky,  Department  of  Mathematics,  Lexington,  k'Y  40506,  USA 
Bergen  Scientific  Centre  IBM,  Allegatcn  36,  N-5007  Bergen,  Norway,  on  leave  from  University  of 
Kentucky,  Department  of  Mat  hematics,  Lexington,  KY  40506,  USA 


1.  Introduction 

In  this  paper  we  focus  on  an  inverse  eigenvalue  problem  for  unitary  Hessenberg  matrices  with  pos- 
itive subdiagonal  elements.  Throughout  this  paper  all  Hessenberg  matrices  are  upper  Hessenberg 
matrices.  This  class  of  matrices  bears  many  similarities  with  the  class  of  Jacobi  matrices,  i.e.  real 
symmetric  tridiagonal  matrices  with  positive  subdiagonal  elements.  Any  matrix  in  either  class  is 
normal  and  has  distinct  eigenvalues.  Both  n  x  n  Jacobi  matrices  and  n  x  n  unitary  Hessenberg 
matrices  with  positive  subdiagonal  elements  can  be  parameterized  by  2n  —  1  real  parameters.  This 
is  obvious  for  Jacobi  matrices;  for  unitary  Hessenberg  matrices  this  parameterization  is  described 
below.  Since  unitary  Hessenberg  matrices  with  positive  subdiagonal  elements  are  determined  by 
O(n)  parameters,  one  can  develop  efficient  algorithms  for  this  class  of  matrices.  These  algorithms 
are  analogous  with  algorithms  for  Jacobi  matrices.  For  example,  the  unitary  QR  algorithm,  de- 
scribed in  [Gr2],  has  many  similarities  with  the  QR  algorithm  for  Jacobi  matrices.  Another  example 
is  provided  by  the  divide-and-conquer  methods  that  have  been  developed  for  both  the  tridiagonal 
and  unitary  eigenproblems  [Cu],  [DS],  [GRl],  [GR2]. 

In  this  paper  we  show  another  analogy  of  unitary  Hessenberg  matrices  with  Jacobi  matrices,  namely, 
that  a  unitary  Hessenberg  matrix  H  with  positive  subdiagonal  elements  is  uniquely  determined 
by  its  eigenvalues  and  the  eigenvalues  of  a  unitary  rank-one  perturbation  of  H .  The  matrix  H 
can  be  constructed  using  0(n2)  arithmetic  operations  using  an  inverse  unitary  QR  algorithm. 
Similar  results  for  Jacobi  matrices  are  well-established  [BGl],  [BG2],  [GH].  The  inverse  unitary  QR 
algorithm  is  analogous  with  the  algorithm  described  in  [GH]  for  Jacobi  matrices. 

2.  Unitary  Hessenberg  Matrices  and  Szego  Polynomials 

We  refer  to  a  finite  Schur  parameter  sequence  of  length  n  as  a  sequence  of  complex  numbers  {7j}j_i 
with  |7j|  <  1  for  1  <  j  <  n  and  |7n|  =  1.  Also  define  the  complementary  Schur  parameters  {&j}"~l 
by  a  j  :=  y/l  —  \fj\2.  Associated  with  the  finite  Schur  parameter  sequence  {Aj}"=1  is  a  unitary 
Hessenberg  matrix  H  with  positive  subdiagonal  elements 

H  =  H{lu. .  .,7n-i,7n)  :=  Gx{lx)G2{l2) . .  .Gn_1(7n_1)Gn(7n),  (2.1) 

where  the  Givens  reflector  Gj(jjj)  is  the  identity  matrix  of  appropriate  size  except  for  the  2x2 
principal  submatrix 

rj   i  +  i 


Gj 


j     j+l 


-7j      <jj 


with  the  bar  denoting  complex  conjugation.  The  matrix  on  the  right  in  the  product  (2.1)  is  defined 
by  Gnifn)  :=  diag[l,  1, . . .,  1,  -7„].  The  nonzero  entries  of  H  =  [Vj,k]jk=i  are  t^ien  given  by 
rjj+i  j  :=  od  and  n^k  :=  -7j_1ai(7j+1  . .  .Ok-ilk  for  1  <  j  <  k,  where  70  :=  1. 

It  is  easy  to  see  that  every  n  X  n  unitary  Hessenberg  matrix  H  =  [Vj,k]1jtk-i  with  Vj+iJ  >  0  is 
uniquely  determined  by  a  finite  Schur  parameter  sequence  of  length  n.  In  fact,  the  Schur  parameters 
{jj}]=i  and  the  complementary  Schur  parameters  {oj}*-!]  can  be  determined  from  //  by 


a3  =  rjj+ij,        1  <  j  <  n;        jj  =  -r)i,j/°i°2  ■  -'^j-i,        1  <  j  <  n. 

Hence,  we  have  a  one-to-one  correspondence  between  Schur  parameter  sequences  of  length  n  and 
nxn  unitary  Hessenberg  matrices  with  positive  subdiagonal  elements.  This  Schur  parameterization 
of  unitary  Hessenberg  matrices  with  positive  subdiagonal  elements  shows  that  these  matrices  are 
determined  by  2n  -  1  real  parameters:  the  real  and  imaginary  parts  of  jj  for  1  <  j  <  n,  and  the 


argument  of  fn.  To  avoid  numerical  instability,  however,  we  also  retain  the  complementary  Schur 
parameters. 

Let  //  =  H„  :=  #(71,  ■  ■  ■  ,7n-i,7n),  and  let  Hk  :=  Gi (71)^2(72)  •  •  •Gk-i('Tk-i)Gk('l(k)  be  the 
leading  principal  sub-matrix  of  II  =  II n  of  order  k.  Introduce  the  functions 

4>k{\):=eT,(XI-Hk)-le^ 
lMA):=det(A/-tffc), 
jr*(A)  :=  det(AJ  -  H'U  ), 

where  e\  :=  [1,0, ...  ,0]T,  and  H't!_1  denotes  the  trailing  principal  submatrix  of  II k  of  order  k  —  1. 
Thus,  <f>k(X)  —  7rA:(A)/t/'fc(A).  The  following  proposition  can  be  verified  by  induction. 

Proposition  2.1.  The  polynomials  V')fc(A)  and  ^(A),  A;  >  0,  satisfy 


<Ma; 


7Tfe(A) 

**(A) 


A 

7  a- 

^fc-l(A)       TTfc-l(A) 

Vo(A) 

7T0(A 

7fcA 

1 

V'a-i(A)      7Tfc_l(A) 

) 

V'o(A) 

7T0(A 

1/A 

0 


It  follows  that  for  each  A*,  the  polynomials  V'A-(A)  :=  A*V^.(1/A)  are  obtained  by  reversing  and 
conjugating  the  coefficients  of  0/-(A).  From  the  initial  conditions  we  recognize  va-(A)  to  be  the  A'th 
Szego  polynomial  determined  by  the  Schur  parameter  sequence  {7j}j=1-     n 

The  Szego  polynomials  {ipk}k=o  are  orthogonal  with  respect  to  a  discrete  measure  on  the  unit 
circle.  This  measure  assigns  a  positive  weight  ujk  to  each  zero  A^  of  «/'n(A).  These  weights  are  the 
numerators  in  the  partial  fraction  decomposition 


<MA)=:]T 


uk 


k=\ 


A-  X, 


see  [Grl]  for  details. 

Let  H  =  UAU*,  with  A  =  diag[A1?  A2, . . . ,  An]  and  U  unitary,  be  the  spectral  resolution  of 
//  =  H(fi, . . .  ,7n_i,  7n),  where  *  denotes  transposition  and  complex  conjugation.  Let  u  :—  I'1 ej 
be  the  vector  containing  the  first  components  of  the  eigenvectors  of  H .  Since  II  has  nonzero  sub- 
diagonal  elements,  every  entry  of  u  =  [v\,  t>2, . . . ,  iJn]T  is  nonzero,  and  we  normalize  U  so  that  each 
vk  >  0.  Then 

n  o 

0n(A)=   U*{XI-\)-lU 


n 
k=l 


A, 


Hence,  the  weights  u>k  —  v\  are  guaranteed  to  be  positive,  and  ^£=1  w*  —  1- 

3.  The  Inverse  Unitary  QR  Algorithm 

Given  n  distinct  unimodular  complex  numbers  {A/j£=1  and  associated  positive  weights  {u£}jj_j, 
we  can  construct  a  unitary  Hessenberg  matrix  //  with  the  Ajt  and  Vk  equal  to  the  eigenvalues  and 
first  components  of  the  corresponding  eigenvectors  of//,  respectively.  This  construction  is  achieved 
using  an  inverse  QR  algorithm,  which  is  analogous  with  the  procedure  of  [Gil]  for  real  symmetric 
tridiagonal  matrices. 

The  required  Hessenberg  matrix  is  obtained  by  performing  a  sequence  of  elementary  unitary  simi- 


laritv  transformations  to  transform  the  matrix 


A 


to  a  Hessenberg  matrix 


11 


without 


using  or  changing  the  arbitrary  entry  8.  Then  //  =  H(ji, . .  .,7n-i,7n)  has  the  desired  eigenvalues 
and  associated  eigenvectors. 

The  idea  is  to  build  the  Hessenberg  matrix  by  adding  weight-abscissa  pairs  one  at  a  time.  Suppose 
that  we  have  constructed  the  unitary  Hessenberg  matrix  Hm  :=  #(7i, . .  •  j7*n-i>7m)s  f°r  some 
m  <  n,  corresponding  with  the  weight-abscissa  pairs  {(u;^,  Afc)}™=1.  Let  gq  :—  (Y^k=\  ^k)1^2  and 
assume  that  the  first  components  of  the  eigenvectors  of  Hm  are  {(t>fc/co)}™=1-  In  order  to  add  the 
weight-abscissa  pair  (v2,  A)  and  construct  the  corresponding  [m  +  1)  X  (m  +  1)  unitary  Hessenberg 
matrix  H'm+l  :=  H(X'1,  A2,  ■  ■  . ,  X'm+\),  we  perform  a  sequence  of  unitary  similarity  transformations 
to  put  the  (777.  +  2)  x  (777  +  2)  matrix 


AM:= 


6 

V 

aQe\ 

V 

A 
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0 

Hm 
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is  a  Hessenberg  matrix  with  a  trailing  principal  (m  +  1)  x  (m  +  1)  submatrix,  which  is  both  unitary 
and  of  Hessenberg  form.  On  the  completion  of  the  similarity  transformation  of  H^l\  we  obtain 


HWG*2{aQ)  = 


(3.1) 


The  circled  element  in  (3.1)  forms  a  "bulge",  which  is  to  be  chased  down  along  the  subdiagonal  in 
order  to  obtain  a  matrix  of  Hessenberg  form.  Define  G^ai)  so  that  G^a^H^Gofao)  =:  H^  is 
a  Hessenberg  matrix.  Then  H^G^(qi)  has  a  bulge,  which  we  annihilate  by  multiplying  from  the 
left  with  6*4(0:2).  Proceeding  in  this  manner,  we  ultimately  chase  the  bulge  off  the  bottom  of  the 
matrix,  and  obtain  the  Hessenberg  matrix  ^^m-1;G^_1(7m_3),  which  is  unitarily  similar  to  H^\ 
The  trailing  principal  (m  +  1)  x  (777  +  1)  submatrix  of  iT^Tn-1)G^_1(7m_3)  is  unitarily  similar  to 
a  unitary  Hessenberg  matrix  with  positive  subdiagonal  elements.  The  latter  matrix  is  the  desired 
Hessenberg  matrix  H'm+l  =  ff(7i,72>  •  •  •  >7m+i)- 

This  procedure  for  adding  a  weight-abscissa  pair  to  Hm,  if  implemented  by  directly  manipulating  the 
elements  of  the  matrices  H^k\  for  1  <  k  <  m,  would  require  0(m2)  arithmetic  operations.  However, 


we  note  that  for  each  k  the  trailing  (m-j-1  )x(m+l)  principal  submatrix  of  I~I(h>  is  unitary  and  of  Hes- 
senberg  form,  and,  therefore,  is  unitarily  similar  to  a  unitary  Hessenberg  matrix,  denoted  by  #„j+ii 
with  positive  subdiagonal  elements.  Hence,  we  can  carry  out  the  similarity  transformations  by  ma- 
nipulating the  Schur  parameters  of  the  matrices  IIm+{  -  #(7{,72>  •  •  ->7a;5  ^kak,~fk+i,  •  •  -,7m)- 
This  gives  rise  to  a  method  that  requires  only  O(m)  arithmetic  operations  in  order  to  add  a  weight- 
abscissa  pair  to  Hm.  We  refer  to  this  method  as  the  inverse  unitary  QR  algorithm  because  of  its 
relationship  with  the  unitary  QR  algorithm  presented  in  [Gr2]. 

Inverse  Unitary  QR  Algorithm:  adding  a  weight-abscissa  pair  (r2,  A),  where  |A|  =  1.  v  >  0. 


<r, 


o  •-  \/ao  +  v2\ 
a0  :=  -v/a'oi  (30  :=  aQ/(r'0; 
for  k  :=  1,2, . . . ,  m 


o\.  :-  fik-iJal  4-  |Qfc_i  +  7fcAfc-2afc_i|2; 


k 

Pk 

L7l- 
7m+i  :=  A7m; 


=  fa-iH^-i  +  7fcA*   2aA._1)/a[.; 
=  0k-i<rk/<r'ki 


=  /?Li7fc  "  Afc"2a 


A--1  i 


Thus,  the  inverse  unitary  QR  algorithm  can  be  used  to  construct  the  unitary  n  x  n  Hessenberg 
matrix  H  =  #(71, . .  -,7n-i?7n)  from  its  eigenvalues  and  the  first  components  of  its  normalized 
eigenvectors  in  0(n2)  arithmetic  operations. 


4.   An  Inverse  Spectral  Problem 


E 


u'k 


Let  H  =  #(7i,---,7n-i,7n)  and  4>n(\)  =     n         - 

v'n(A)       ^  A  -  AA- 

and  consider  the  polynomials 

XkW  ■=  (1  -  a)A7r,-(A)  +  o^-(A),      0  <  k  <  n. 


Let  a  :=  elT  for  some  0  <  r  <  27r. 


(4.1) 


Proposition  4.1.    The  polynomials  {\'a-(A)}£_0  are  the  monk  Szego  polynomials  corresponding 
with  the  Schur  sequence  {o7a.}£=i- 

Proof.  The  definition  (4.1)  and  Proposition  2.1  yield 
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e*(A) 


A       7fc 
7JfcA      1 


\a-i(A)1 

[\o(A)l 

r  1 

e*-i(A)j 
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Leo(A)  j 

a 

where  ^-(A)  :=  (1  -  a)A7ffc(A)  +  a^fc(A).  On  setting  X>(A)  =  a^.(A),  we  obtain 


Xfc(A) 

A 

Otlk 

Xfc-i(A) 

\o(A 

U(A) 

Q7*  A 

1 

U-i(A) 

? 

[\o(A 

An  immediate  consequence  of  Proposition  4.1  is  that  each  zero  of  \„(A)  lias  unit  modulus.  When 
a  =  -1  the  polynomials  {\jt(A)}£_0  are  known  as  the  Szego  polynomials  of  the  second  kind 
corresponding  with  {7fc}£=i> 


Proposition  4.2.  The  zeros  {fJ-k}k-o  of \n(A)  strictly  interlace  the  zeros  {A^.}^_0  of  ipn(X)  on  the 
unit  circle.  Moreover,  Ylk=i(^k/^k)  —  Q- 

Proof.  Let  A  :=  et9  for  some  0  <  6  <  2x,  and  let  a  =:  etT  for  some  0  <  r  <  2n.  Then 
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Xn(A)         V"^         A  -  qA* 


MA):=^  =  E 


uk- 


V'n(A)       ^        A  -  Xk 
A      1  -  AAfc  +  a(l  -  AAfc) 


fc=i 


|A-A,| 


_^9    .>/2  ,  Re(c-/2(1-AA,)) 


=  2>,T/2w* 


fc=l 


|A-A,|2 


Thus,  Im(e-IT/2$n(e^))  =  0  for  0  <  6>  <  2?r.  Let  Afc  =:  e*'** .  We  may  assume  that  0  <  0X  <  62  < 
■  ■  ■  <  9n  <  2tt.  Then 

cos(r/2)  -  cos(r/2  +  0*  -  #) 


c  *n(£      )  =    >      iWjt ; ^ o 

&        fc(l_cos(0, -0))2+sin2(#, -0) 


fc=l 

cos(r/2)  -  cos(r/2 

UJl. ■ 

1  —  cos(^/c  — 


cos(r/2)  -  cos(r/2)  cos(6k  -  6)  +  sin(r/2)  sin(0fc  -  0) 

ujk 


=  ]Tu?fc 


=  ^  uk  (cos(r/2)  +  sin(r/2)cot((0fc  -  0)/2)) . 
From  0  <  r  <  27r  and  uk  >  0  for  all  fc,  it  follows  that 

e-!r/2  — $n(e^)=  -sin(r/2)^^,/sin2((^-^)/2)>0     for     0^  9kt  1  <  k  <  n. 

Thus,  $„(eie)  -+  -oo  as  0  \  Ok,  and  $n(el(?)  -+  oo  as  0  /  8k+1.  This  shows  that  9  -►  $n(eie)  has 
precisely  one  zero  in  ]6k,0k+\[.  Consequently,  the  zeros  fij  of  $n(A)  strictly  interlace  the  Xj  on  the 
unit  circle.  The  second  statement  of  the  Proposition  follows  from  the  fact  that  Xn(0)/V'n(0)  —  <*■     D 

Let  Xk  =:  ei6k  and  fik  =:  e*"*  for  0  <  dk,vk  <  2tt.  Then  we  have  ELiK  ~  ^*)  =  r-  We  ma>' 
assume  that  the  arguments  have  been  ordered  so  that 

0  <  01   <  VX  <  #2   <  •  •  •  <  "n-X    <  On   <  »n   <  #1  +  2tT.  (4.2) 

Proposition  4.3.  With  the  above  notation,  and  the  ordering  (4.2),  the  weights  u>k  are  given  by 

ft  ain(to  -  ^-)/2) 
1         j=i 

U!k   = . 

sin^2)  ft  -b(»  -  »»)/2) 

Proof.  We  have 


m  mm       n(A*-w) 

1  -  QjAfcWjt  =    lim    A  -  Afc)  =  =  — 

^n(A)       ^(A,)         n(Afc_Ai) 


since  Xn(^)  an<^  ^n(^)  are  monic.  Thus. 

f[(ei6k  -eivi)        n(l-e*<"i-^))  ft  sin((^  -  6>fc)/2) 

(l-ap,  =  e-^^-  -  =  ^-  -  =  -2^if  _,      (4.3) 

Yl  (ctff*  _  cWy)         ft  (1  -  c*'(*i-«*))  ft  sin((^  -  0fc)/2) 

because 

1  -etp  =  -2isin(/5/2),  (4.4) 

and  $^^=1(i/j  —  #,)  =  r.  The  formula  for  the  weights  now  follows  by  substituting  (4.4)  with  /3  :=  r 
into  (4.3).     □ 

We  can  now  state  the  inverse  spectral  problem  and  its  solution. 

Problem:  Given  two  sets  of  n  mutually  interlacing  points  on  the  unit  circle  {A^}^=1  and  {/u-}"=1. 
determine  the  unique  unitary  Hessenberg  matrix  H  =  //(71, . . . ,7n-i)7n)  and  |o|  =  1  such  that 
the  spectrum  of//  is  {A/.}£=1  and  the  spectrum  of  II(a~f\, . . . ,  a7n_i,  Q7„)  is  {/U-}£=i- 

Solution:  Let  a  :  —  Yl^ii Vk/^k),  calculate  the  weights  by  Proposition  4.3.  and  use  the  inverse 
unitary  QR  algorithm  to  construct  H  —  //(71, . . . ,7n-i,7n)- 

References 

[BGl]  D.  Boley  and  G.H.  Golub,  Inverse  Eigenvalue  Problems  for  Band  Matrices,  in  Lecture  Notes 
in  Mathematics  630,  Numerical  Analysis,  Proceedings  of  the  Biennial  Conference  Held  at 
Dundee,  1977,  G.A.  Watson,  ed.,  Springer- Verlag,  New  York.  1978,  pp.  23-31. 

[BG2]  D.  Boley  and  G.H.  Golub,  A  Survey  of  Matrix  Inverse  Eigenvalue  Problems,  Inverse  Problems 
3,  (1987)  595-622. 

[Cu]  J.J.M.  Cuppen,  A  Divide  and  Conquer  Method  for  the  Symmetric  Tridiagonal  Eigenproblem, 
Numer.  Math.  36,  (1981)  177-195. 

[DS]  J.J.  Dongarra  and  D.C.  Sorensen,  A  Fully  Parallel  Algorithm  for  the  Symmetric  Eigenproblt  m, 
SIAM  J.  Sci.  Stat.  Comput.  8,  (1987)  sl39-sl54. 

[Gil]  W.B.  Gragg,  Positive  Definite  Toeplitz  Matrices,  the  Arnoldi  Process  for  Isometric  Operators, 
and  Gaussian  Quadrature  on  the  Unit  Circle  (in  Russian),  in  Numerical  Methods  in  Linear 
Algebra,  E.S.  Nikolaev,  ed.,  Moscow  University  Press,  Moscow,  1982.  pp.   16-32. 

[Gr2]  W.B.  Gragg,  The  QR  Algorithm  for  Unitary  Hessenberg  Matrices,  J.  Comput.  Appl.  Math. 
16, (1986)  1-8. 

[GH]  W.B.  Gragg  and  W.J.  Harrod,  The  Numerically  Stable  Reconstruction  of  Jacobi  Matrices 
from  Spectral  Data,  Numer.  Math.  44,  (1984)  317-355. 

[GR1]  W.B.  Gragg  and  L.  Reichel,  A  Divide  and  Conquer  Algorithm  for  the  Unitary  Eigenproblem. 
in  Hypercube  Multiprocessors  1987,  M.T.  Heath,  ed.,  SIAM.  Philadelphia,  1987,  pp.  639-647. 

[GR2]  W.B.  Gragg  and  L.  Reichel,  .4  Divide  and  Conquer  Method  for  the  Unitary  and  Orthogonal 
Eiejenproblems,  BSC  Report  88/31.  Bergen  Scientific  Centre,  Bergen,  1988. 


DISTRIBUTION  LIST 


DIRECTOR  (2) 

DEFENSE  TECH.  INFORMATION 
CENTER,  CAMERON  STATION 
ALEXANDRIA,  VA   22314 

DIRECTOR  OF  RESEARCH  ADMIN. 
CODE  012 

NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA   9  3  943 

LIBRARY  (2) 

CODE  0142 

NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA   93  94  3 

DEPARTMENT  OF  MATHEMATICS 
CODE  53 

NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA   9  394  3 

CENTER  FOR  NAVAL  ANALYSES 
44  01  FORD  AVENUE 
ALEXANDRIA,  VA   22302-0268 

PROFESSOR  WILLIAM  GRAGG  (15) 

CODE  53Gr 

DEPARTMENT  OF  MATHEMATICS 

NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY,  CA   9  3  943 

NATIONAL  SCIENCE  FOUNDATION 
WASHINGTON,  D.C.   20550 


DUDLEY  KNOX  LIBRARY 

3  2768  00354737  3 


